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STABILITY FOR QUANTITATIVE PHOTOACOUSTIC 
TOMOGRAPHY WITH WELL CHOSEN ILLUMINATIONS 


GIOVANNI ALESSANDRINI, MICHELE DI CRISTO, ELISA FRANCINI, 
AND SERGIO VESSELLA 

Abstract. We treat the stability issue for the three dimensional inverse imag¬ 
ing modality called Quantitative Photoacoustic Tomography. We provide uni¬ 
versal choices of the illuminations which enable to recover, in a Holder stable 
fashion, the diffusion and absorption coefficients from the interior pressure 
data. With such choices of illuminations we do not need the nondegeneracy 
conditions commonly used in previous studies, which are difficult to be verified 
a-priori. 


1. Introduction 


In recent years there has been an increasing interest towards imaging methods 
which combine different physical modalities of interrogation which are known as 
hybrid, or coupled physics, inverse problems. 

In this paper we concentrate on Photoacoustic Tomography (PAT) that couples 
the (high contrast) optical tomography with the (high resolution) ultrasound waves. 
The first step of this imaging technique consists of reading the boundary response 
to an acoustic signal, in order to reconstruct the absorbed energy distribution in¬ 
side the biological tissue under inspection. We assume this first part as already 
performed and we focus on the second step of the procedure, which consists of 
reconstructing the absorption and diffusion coefficients from measurements of the 
absorbed energy distribution obtained in the previous step. We refer to [RGZ13] 
for an extended bibliographical review on this problem known as Quantitative Pho¬ 
toacoustic Tomography (qPAT). 

Let us denote by C K" the body enclosing the biological tissue under inspec¬ 
tion. Note that the physically significant space dimension is n = 3, but, since we 
shall need also to refer to model cases when n = 2, we shall leave n > 2 undeter¬ 
mined. We denote by u(x) the photon density at the point x £ fl. Then u(x ) solves 
the following boundary value problem 

DX7u) + cru = 0 in Q, 



( 1 . 1 ) 


where D = D{x) > 0, a = cr(x) are the diffusion and absorption coefficients, 
respectively, and g is the illumination source prescribed on the boundary. 

The goal of qPAT is to recover information on the coefficients D , <j from the 
knowledge of the absorbed energy 


H(x) = cr(x)u(x), x £ fl, 


possibly repeating the experiment with different profiles of the illumination g. 
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Note, incidentally, that a more accurate model would require the introduction 
of the additional multiplicative unknown parameter r(x) called the Griineisen co¬ 
efficient. Here we adopt, for simplicity, the commonly used convention of assuming 
r = 1. See iBRlll for a discussion on this issue. 

This problem has been considered in [BUlOj . where a uniqueness result with 
two measurements is proven. The authors assume D G C k+2 , a G C k+1 with 
k > 1. They also present a Lipsclritz stability theorem. This result has been later 
improved in [BRlTl also providing a numerical reconstruction procedure. More 
recently, [ RGZ13I have considered the same problem when the prescribed illumi¬ 
nation is modeled by a Robin boundary condition, rather than a Dirichlet one. A 
reconstruction method for the full nonlinear inversion is also treated. 

All the above quoted results rely on a nondegeneracy condition which can be 
illustrated as follows. If iq, U 2 are two solutions to mi) corresponding to two 
different illuminations < 7 i, <72 > 0, then it is a quite well-known fact that the quotient 

U-2 

u = — 

Ul 

satisfies an elliptic equation in pure divergence form 

f div(aVu) =0 in H, 

[ u = g on dfl, 

where g is the ratio ^ (see Proposition 14.II below - ). It is also easy to see that the 
solution of the qPAT problem boils down to solve the inverse problem of finding a 
given u and the boundary values a\ ga . This is a relatively easy task if g is chosen 
in such a way that the following nondegeneracy condition holds: 

(1.2) |Vm| > 0, everywhere in Q. 

When n = 2 there exists a well-established criterion which enables to choose the 
Dirichlet data g (independently of a!) so that (11.21) holds ( IA861 . [AM94] - ). Such 
a criterion is unimodality. That is, roughly speaking that the graph of g has one 
single peak of maximum points, one of minimum and is monotone in between. 

In dimension n > 3 complex valued solutions satisfying m can be constructed 
by the method of Complex Geometrical Optics [BUIOI , but their boundary data do 
depend on the interior values of the (unknown) coefficient a and thus they cannot 
be a-priori chosen. 

Real valued solutions which locally satisfy |Vu| > 0 can also be constructed (see 
|GW75i . ( BU131 Theorem 4.7]), but still they depend on the unknown coefficient a. 

Indeed, there are reasons to believe that, when n > 3, it does not exists any 
Dirichlet data g such that is satisfied for every a. See, for related discussions, 
[AN15I . [CII5]. 

The principal aim of the present paper is to treat stability even when the above 
stated nondegeneracy condition may be violated. 

We shall show that stability can be obtained with essentially arbitrary illumi¬ 
nations gi 7 <72 imposing only one constraint on their ratio g = Namely that g 
satisfies a condition of unimodality adapted to the (n — l)-dimensional boundary 
dfl. This condition shall be made more precise in Definition 12.31 

Our strategy shall be as follows. In dimension n > 3, for a fixed g, we cannot 
assure the nonvanishing of |Vu| throughout B, but, under reasonable assumptions, 
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it is possible to keep under control the vanishing rate in the interior. This will be 
the content of Lemma 13.31 

On the other hand, assuming unimodality we can make sure that |Vu| > 0 on a 
small neighborhood of <9f l (Lemma 13.41) . 

Next we adapt from j A.86] a weighted stability estimate on the coefficients a in 
terms of u and of a\ gn . 

Using the previously mentioned estimates on the vanishing rate of |Vu| and 
a suitable interpolation inequality [Si 151 . we arrive at an (unweighted) stability 
estimate of Holder type for a (Theorem 13.51) . The deduction of stability bounds for 
D and a follows the track of well-known arguments, see for instance |RGZ13[ . 

Let us emphasize that most of the present effort is devoted to two main goals: 

(1) To avoid the nondegeneracy condition m- 

(2) To make precise (but feasible) a-priori assumptions which guarantee a quan¬ 
titative, concrete, evaluation in our stability estimates. 

It is our belief that the present approach can be useful also to other, more complex, 
hybrid inverse problems, where analogous issues of nondegeneracy arise. 

The paper is organized as follows. In the next Section [2] we provide the main 
assumptions and we state our main result (Theorem 12.4[) . The proof of it is based 
on some auxiliary propositions, given in the subsequent Section [3j and is presented 
in Section [4] 


2. Assumptions and Main Result 

We assume to be a C 2 - smooth, bounded domain in R ra diffeomorphic to 
the unit ball Hi(0). More precisely, from a quantitative point of view, we assume 
that there exists a diffeomorphism F of class C 2 such that, for given constants Q o, 


Q i > 0, 


(2.1a) 

F : Ri(0) GA n, 

(2.1b) 

II- f IIc 2 (b 1 (o)) < Q o, 

(2.1c) 

I F(x) - F(y )| > \x - y\, for every x,y € Ri(0) 

Wi 


The constants Q o, Q i, shall be part of the a-priori information that shall be used 
in our quantitative stability estimates. 

We are interested in recovering the unknown parameters D and a, by performing 
two measurements, that is prescribing two data tq and <?2 on dtt and measuring 
the corresponding internal pressure fields. In particular we establish a continuous 
dependence of the unknown parameters from the measured data. 

Given constants Ao, Ai, Eq, E±, qo, qi, q 2 > 0 (which also shall be part of the 
a-priori information), we consider unknown coefficients D and a such that 


(2.2) 

D G W ll00 (f2), 

a G W 1,00 (fl), 

and 



(2.3) 

Aq 1 <D< Ao, 

for every igU, 

(2.4) 

A^ 1 < o < Ai, 

for every iGd, 
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(2.5) ||I?||vK 1 >“(n) < IkllvK^^cn) < E\. 
The boundary data we choose are functions g\ and g 2 such that 

(2.6) gi £ C 2 {dn), WgiWcHdQ) < £to for * = 1,2, 


(2.7) Mi 1 < 9i{ x ) < Mi, for every x £ fl, * = 1,2. 

Moreover, denoting by 

32 , - 1 /' 

g = — and g = j—-y / 3 , 

3i Jon 

we assume that 


( 2 - 8 ) II 5 - 3 IU 2 (dfi) > M2 1 - 

Remark 2.1. Let us emphasize that (12.811 represents constructively the assumption 
that the illuminations g\, g 2 are linearly independent. 


Remark 2.2. Note that, if gi, g 2 satisfy assumptions m and (12.81) . then the so- 
called frequency function associated to g 


(2.9) 


F[9] := 


II3-3II 


if!/2(on) 


Il3 - 3lU 2 (an) 

is bounded by a constant depending only on go, g 2 and Q\. 


We shall see that it is convenient to assume that the ratio q = — of the illu- 

a gi 

minations has a specific behaviour which is expressed in the following definition. 
This is a form of monotonicity assumption, which needs however to be specified in 
quantitative fashion. 


Definition 2.3. Given m, M, 0 < m < M and a continuous, strictly increasing 
function u 1 : R + —>• R + , such that w(0) = 0, we say that a function g £ C ,1 (9f2,R) 
is quantitatively unimodal if 

(2.10) m < g(x) < M, for every x £ <9f l, 
the subsets of dfl 

(2.11) T m = {x £ dfl : g(x) = to} and Tm = {x £ dfl : g(x) = M} 

are connected and non-empty, possibly reduced to single points, and, for every x £ 
dfl \ (T m U Tm) such that dist(x, T m U Tm) > S, we have 

(2.12) |V Tff (x)| >u{6), 
where Vt denotes the tangential gradient. 


Also m, M , and ui shall be part of the a-priori information. 

In our stability result, we compare two different sets of diffusion and absorption 
coefficients, that we denote by D^\ <r f1 ^ and D^ 2 \ . 

Let m-A for i, j = 1, 2, be solution to 


(2.13) 


-div (i)(j)Vu| i) ) + cr( j '>u < ? ) = 0 in O. 

= gi on dCl. 




STABILITY FOR QPAT 


5 


(i) 

We emphasize that, in the notation uf, the superscript is associated to the un¬ 
known parameters D^\ whereas the subscript is associated to the illumination 
9i- 

The available measurements which represent the internal pressure fields gener¬ 
ated by the absorptions of photons energy are given by 

(2.14) H\ j) =a < ~ j) u ( f ) . 

We prove the following (uniqueness and) stability result: 


hP - hP 


Theorem 2.4. Let all the assumptions stated above be satisfied. If 

(2.15) 

and 

(2.16) 


l 2 ( n) 
D (1) -Z? (2) 


then we have 
(2.17) 


£)(1) _ £)(2) 

+ 

cr^ 1 ' — a 1 ' 2 ' 1 


L 2 (Q) 



< e, fori = 1,2, 

<e', 

< C (e + ep , 

l 2 ( n) 


L°°(d Q) 


where C and 6 G (0,1) only depend on the a-priori information Qo, Qi, Ao, Ai, 
Eq, Ei, go, hi, p. 2 , in, M and to. 


3. Auxiliary results 

The proof of Theorem 12.41 is based on a result concerning stable reconstruction 
of the main coefficient of a second order elliptic equation from the knowledge of 
internal values of one of its non constant solutions. 

Theorem 3.1. Let Li be diffeomorphic to the unit ball. Let a and b G W 1 , 00 (f2) 
such that 

(3.1) Cp < a(x), b(x) < Co, for every x G Cl, 
and 

(3.2) |Va(x)|, |V 6 (x)| < C\, for almost every x G Cl. 

Let g and k G C 2 (dCl), with 

(3-3) ||ff||c 2 (dn), ||fc||c 2 (on) < C 2 . 

Assume g satisfies assumption (12.81) and 

(3.4) g(x) > Cp for every x G LI. 

be the unique solutions to boundary value problems 

div (a(x)Vu(x)) =0 in Cl, 
u = g on dCl, 

div (b(x)Vv(x)) = 0 in Cl, 
v = k on dCl. 

Given Cl' CC Cl with dist{dCl , Cl') > do and 9 G (0,1/2) there are positive constants 
C and /3, depending only on Co, C\, C 2 , H 2 , C 3 , Qo, Q 1 , do and 9, such that 

a - &||z~(n') < & (jl u _ v \\ e L 2 (n) + ll a — &IU°°(9n)) • 


Let u and v in W 1 ’ 2 (Cl) 

(3.5) 
and 

(3.6) 


(3.7) 
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In order to prove Theorem l3.il we need the following lemmas. For d > 0 we will 
denote by fid = {x £ : dist{x,d£i) > d}. 

Lemma 3.2. With the same assumptions of Theorem \3.H for every p > 0 and for 
every x £ fl^ p 

[ |v M | 2 > c p [ |v«| 2 , 

J B p {x) JQ 

where C p depends on Co, C\, Qo, Qi, F[g] and p only. 

Proof. This Lemma corresponds to Theorem 4.2 in |AMR03j for solutions of Dirich- 
let boundary problem instead of Neumann boundary problem. The proof follows 
the same path by taking uq = u — u where u = |fl| _1 Jq u - The only difference sits 
in getting from estimate (4.28) in [AMR03j to an estimate in terms of the frequency 
function F\g). In this case, by standard elliptic estimates, 

l|Vuo||i 2 (n) = ||V(it - g)\\L 2 (n) < C\\g - g\\ H l / 2 (an)^ 
where C depends on Q 0 , Q\ and Cq, and 

||«o|U 2 (9n) = II<7 — w||i 2 (dn) > \\g - 5|U 2 (an), 
hence, from estimate (4.28) in fAMR03I . we get 

l|Vuo|| L 2( n ) ^ ll VM o|||2 ( n) ^ / llff - g\\m/^dn) \ 2 _ rF , -p 

ll«oll^ (n) " IKH£ W - ' V h-ahnon) J m ' 

The rest of the proof is as in [AMR03] . □ 


Lemma 3.3. With the same assumptions of Theorem Iff. 1 1 given ii' CC S2 with 
dist(dfl,Tt') > do, there exist positive constants K\ and K 2 > 1 and d\ < do, 
depending only on Co, C\, C 2 , Pi, C 3 and do, such that, for every xq £ ft' and 
r < di, 


(3.8) 


/ |Vu|2 - IT' 

lB r (x 0 ) F 2 


Proof. As in the proof of Theorem 4.3 in [AMR03I . starting from the doubling 
inequality by Garofalo and Lin ( |GL 86 j 1 and using Caccioppoli and Poincare in¬ 
equalities, we can show that there is a constant 2 di < do, depending only on Co 
and Ci, such that we get, for r < p < d\ 

(3.9) f |Vu | 2 < C [ |Vu| 2 , 

JB p (x 0 ) \ 1 J J B r (x0) 


where C and K depend on Co and Ci, and I\ depends also, increasingly, on 


(3.10) 


JV(di) = 


d U B 2 dl tx o) \^ u \ 2 

fs 2dl (x 0 )( U ~ Ur ) 2 


where u r = | B J g ^ x ^ u. Now, in order to estimate N(di) from above in terms 
of a-priori information, we can use Caccioppoli inequality to get 


N(di) < 


C c fj 


B 2 d 1 (xo) 


|Vu| : 


4 l 


B dl (xo) 


|Vu| 


(3.11) 
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(3.12) 


By (13.111) and Lemma 13.21 we have 

Cc In l^l 2 Cc 


N(di) < 


4 C dl /JV U | 2 4<V 

By (13.91) with p = d\, (|3.12l) and Lemma 13.21 again, we finally get 

Ki t 


(3.13) / |Vu | 2 > C _1 . 

JB r {x o) V«i 


r 2 ^ |V«| 2 >^ / |Vm| 2 . 


J Bi ± (x o) 

Now, by assumption (12.81) , trace estimates and Poincare inequality we get 
(3.14) p 2 2 < ||g — g\\ 2 H i /2 < C\\u — 5 ||^i(o) < C1I Vu||| 2 (n), 

hence, by combining (13.131) and (13.141) . we finally get (13.81) . 


□ 


Proof of Theorem \3.1l The first step in the proof relies on Lemma 2.1 contained 
in [A 86 j . Although such Lemma 2.1 is stated in a two-dimensional setting, its 
validity can be extended in a straightforward manner to any dimension. 

By Lemma 2.1 of j A~ 86 ] . for any 9 £ (0,1/2) there is a constant K 0 > 0, depending 
only on Co, Ci, C 2 , Q 0 , Q 1 and 0 , such that 

(3.15) J r | a - b\ |Vw| 2 < K 0 (||u - v\\ e L 2 (n) + ||a - &|U°°( 9 fi)) . 

We now reproduce an argument due to Sincich |Sil51 Proposition 4.9]. Let us 
set (f> = a — b and let Xg £ 12 ' such that 


(3.16) 


\(/>{xo)\ = max|<£(x)|. 
O' 


Since a and b satisfy assumption (E2D, 

(3.17) |</>(a7o)I < |0(aO| + 2Cir, for every x £ B r (x 0 ), with 0 < r < dg. 

Multiplying (13.171) by |Vu(x )| 2 and integrating with respect to x on B r (x 0 ), we 
get 


(3.18) \4>{x 0 )\ 
hence 

(3.19) 


/ \Vu(x)\ 2 dx < 

' B r (x 0 ) 


' B r (x 0 ) 


(x)\\Vu(x)\ 2 dx + C\r 


/ |Vw(x)| 2 dx, 

/ B r (x 0 ) 


ljt , „ / f Br (xo) \ 0 ( x )\\V u (x)\ 2 dx 

I</>(x 0 )| <- , -h2Cir. 


f Br {x 0) l Vu (*)l 2 dx 

By (]3.15D and (13.81) . and by (13.161) we have 

(3.20) max |a(x) — b(x)\ < K 0 K 2 r~ Kl (j|u - w||i 2(n) + ||a - 6 ||z,oo(g n) j + 2Cir. 
By choosing an appropriate r £ (0,dg) we get (13.71) . □ 


Let us now show that, by choosing a boundary condition with some additional 
features, we can bound form below the norm of Vu in a neighborhood of the 
boundary. The following Lemma is a variation on themes treated in |AN151 Lemma 
2.8, Theorem 4.1], 
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Lemma 3.4. Let Q be diffeomorphic to the unit ball and let g be quantitatively 
unimodal, according to Definition \2.3[ If u is the unique solution of problem (ESD 
for a coefficient a satisfying assumption & then 

(3.21) Vu| > C, for every x £ f2, with dist(x, dfl) < p, 

where p > 0 and C > 0 depend only on Co, Qo, Q±, m, M and ui. 


Proof. The diffeomorphism F in assumption (12.11) transforms the elliptic equation 
in (13.51) into a similar elliptic equation in £?i(0) with IT 1,00 main coefficient. The 
constant of ellipticity and all the constant appearing in the assumptions ( 12 . 101 ) . 
( 12 . 111 ) and ( 12 . 121 ) shall be changed in a controlled manner only depending on the 
a-priori information. For this reason we assume, without loss of generality, that 

ll = Bi(0). _ 

By regularity estimates for solutions of elliptic equations, the i(0)) norm 

of u is bounded in terms of the a-priori information, hence, for x £ l?i( 0 ) and 
dist{x,T m) < p we have 

u(x) — m > M — m — Cg. 


By choosing g small enough, we get 


By Harnack inequality ( [GT981 Theorem 8.20, Corollary 8.21]), 

M — TYl - 

(3.22) u{x) — m > 0^ --- , for every x £ B\- v (0). 

In particular, if we choose y £ T m and x = (1 — g)y , we get 

u{x) - u{y) > C V M 2 m ■ 

By Hopf lemma 1 ICT981 Lemma 3.4]) we have 

|Vu(y)| > k > 0, for every y £ T m . 


Since we can proceed in the same way on Tm and by using again the regularity 
C 1,/3 of u up to the boundary, we have 

(3.23) |Vu(x)| >k — CS for every x £ dtt, with dist(x,T m UTm) < S. 

By choosing 6 = 6 so that C6^ = k/2, by (13.231) and (12.121) . we have 

|Vu| > max{w(5), k/2}, on dfl. 

By using again the C 1,/3 regularity of u up to the boundary, we get (|3.21D . □ 


Theorem 3.5. Let Q, a, b, g and k be as in Theorem \3.1\. Let us also assume 
that g is quantitatively unimodal. Let u and v in W 1,2 (Tl) be the unique solutions 
to boundary value problems (13.51) and (IffiBl) . Given 8 £ (0,1/2) there are positive 
constants C and (3, depending only on Co, C\, Ci, P 2 , Co, Qo, Q i, do, 9, m, M 
and uj such that 

||a - < C ^||u - u||L2 (n) + ||a - &|| L oo( 0 q)) ■ 


( 3 . 24 ) 
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Proof. The proof follows the same steps of the proof of Theorem 13.11 It suffices to 
extend estimate (13.8[) to every point xg in II. This extension is possible by Lemma 
13.41 and by the regularity assumptions on dfl. □ 


4. Proof of Theorem 12.41 

We proceed as in [ROZ13] (and in pBUiOl ) and show that, for a a fixed set of 
coefficients, the ratio of two solutions (corresponding to different boundary values) 
satisfies a partial differential equation. 


Proposition 4.1. For j 

(4.1) 


1 ,2, the function 

Tj(j) 

U 0) = " 2 


II 


U) 


satisfies the equation 

(4.2) - div (a w = 0 in fi, 

where 

(4.3) = ^2 (^P) 2 = DU) (4 j) ) 2 ■ 

Moreover, 

U (j ) = — on dfl. 

9 1 

Proof of Proposition \4-l\ For the sake of simplicity we drop the superscript ( j ). 
It is an easy calculation to check that, since u, solves equation (12.1311 . then 


-div(D«?V(g))=0, 


hence, by (12.141) and (14.31) 


-div|aV|^))=-div(D«;v(^ 


= 0 . 


□ 


Our aim is to apply Theorem 13.51 to the functions U^ 1 ' 1 and U ^ introduced in 
the previous proposition. First of all, let us show that: 


Claim 4.2. If (12.151) holds, then 


(4.4) 


[7 (1) - t/ (2) 


l 2 ( n) 


< Ce, 


where C depends only on Ao, Ai, /.to, Eg and E\. 
Proof of Claim \f.2\ Let us write 


[hP-h 2 


(4.5) 


C7 (1 ) - U {2) = 


( 2 ) 


) h[ 2) + ( 


h[ 2) - H 1 


(1) \ M 2 ) 


h[ 1 ] h[ 2) 
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We need to get a lower estimate for the functions h[^ and h[ 2 \ By (12.1411 and 
m, it is enough to show that i/f 1 is bounded from below in terms of a-priori 
information. 

By (12.4| and (12.71) . we can apply Maximum Principle and get 

(4.6) 0 < Ui( x) < hi, for every x £ Cl. 

By Theorem 8.33 in [ IGT98I . Viti is bounded in terms of a-priori information. 
Since u[ 3 ^ = gi > gf 1 on dfl (by (12.71) 4. there exists a positive constant d , such 
that 

(4.7) u[ j) > in n \ fl d , 

where, we recall, ild = {% € 12 : dist(x, dfl) > d}. By Harnach inequality f )GT98l 
Theorem 8.20, Corollary 8.21]), 

(4.8) C H inf u[ 3) > sup u[ j) > 

^d/2 Q d /2 ^ 

hence, by (12.51) and (12.71) . 

(4.9) inf = min < inf u\ 3 \ inf > > min 

n l^d /2 si\n d j 


[ Mi Mi 1 


l 2 Ch 2 


M3 ■ 


Inequality (14.91) holds for uf ] as well. By (14.91) and (12.41) . 

(4.10) H[ j) = a^u ( p > (Ai^)" 1 - 
Moreover, by (El and El 

(4.11) #P< Ai/xi. 

By (14.51) . (14.101) . (14.101) and (12.151) we finally have 

(4.12) ||t/ (1) -£/ (2) ||L2(fi) <2hi\\hIz- 


□ 


Proof of Theorem \2.4\ (Conclusion). 

We notice that, by (14.31) . (12.31) . (12.71) . (12.51) and by Theorem 8.33 in ) GT98| . 
coefficients and a^ satisfy assumptions El and El- 
Moreover, by (14.31) . (12.71) and (12.161) . 


(4.13) 

Hence, by Theorem 13.51 


- a™ 


L°°(d n) 


< M?e'- 


(4.14) 


||a (1) - a (2) ||z,oo(n) <C{e + e'f 


where C and 9 depend only on the a-priori information. 

We now proceed as in [RGZ13] . By straightforward calculation it is easy to show 
(i) 

that the function u{ solves 



u[ j) 9i 


(4.15) 


—div ( V 


in H, 
on dfl, 
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from which we can get 
(4.16) 


- diT (“ <1,v (;jr- ^i) ) - Hi " ~ H? ' +iw ((*“ - “ <2> ) v 


By Theorem 8.34 in [GT98I 
(4.17) 


v( 1 ) 


Vu{ 2) 

U 2) J 

L°°(Q) 

K ’) 2 




V< 


( 2 ) 


l°°( a.) 


<C, 


l°°( n) 


hence, by (14.141) . by assumption (12.151) and by Corollary 8.7 in [GT981 ■ since -py — 

U 1 

-j 2 j = 0 on dfl, we conclude that 

1 1 


(4.18) 


(i) ( 2 ) 

l “i 


U[ u 
Finally, by definition (12.141) . 

h[ 1] h[ 


<C {£ + £’)'\ 


vk 1,2 (n) 


- a™ = 


4 2) 


O’(l) 7^(2) 

~ , zj{ 1 ) 

-m- + H i 


( 4 - 19) - - m (« ... 

u\ u\ U] 

hence, by (|2.15)) and (|4.18|h 

(4.20) I|cr (1) - cr (2) || L 2 (n) < fx 3 £ + Ai/ii C (e + e') 9 • 

Since 


A 1 ) .A 2 ) 


(4.21) 


£)( 1 ) _ £)( 2 ) = 


7.(1) 


7.(2) 


(“")' (“?’)' 


by (14.141) . (14.91) and (14.181) 

(4.22) \\ D W - D W\\ L2{n) <C(£ + £') e , 

where C and 9 depend only on the a-priori information. 
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